o 
o 
o 

On: 

in 



Oh: 
a: 



U 



Contribution of Inelastic Rescatterring to 
B — > 7T7T, i\"A^ Decays 

P. Zenczykowski* 
Depi 0/ Theoretical Physics, Institute of Nuclear Physics 



> ; Radzikowskiego 152, 31-342 Krakow, Poland 

in 
§y ; February 1, 2008 

O 

o 



Abstract 

We discuss multichannel inelastic rescatterring effects in B decays into 
a pair PP of pseudoscalar mesons (PP = tttt or KK). In agreement with 
short-distance models it is assumed that initially B meson decays domi- 
k>( \ nantly into jet-like states composed of two flying-apart low-mass resonances 



M1M2 which rescatter into PP. Since from all S-matrix elements (i\S\PP) 
involving PP only some (i = MxM-z) contribute to the final state rescat- 
terring, the latter is treated as a correction only. The rescatterring of res- 
onance pair M1M2 into the final PP state is assumed to proceed through 
Regge exchange. Although effects due to a single intermediate state M\Mi 
are small, it is shown that the combined effect of all such states should be 
large. In particular, amplitudes of B decays into KK become significantly 
larger than those estimated through short-distance penguin diagrams, to 
the point of being comparable to the B — > tttt amplitudes. 
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1 Introduction 

Studies of CP-violation in B decays must involve final state interaction (FSI) ef- 
fects. Unfortunately, a reliable estimate of such effects is very hard to achieve. In 
the analyses of B —* PP decays (P - pseudoscalar meson) only some intermedi- 
ate states, believed to provide nonnegligible contributions, are usually taken into 
account. Many authors restrict their studies to elastic rescatterring P1P2 —>■ P\P<i 
only. In Regge language this is described in terms of a Pomeron exchange. Al- 
though in PP — ► PP quasi-elastic rescatterring at s = m 2 B contributions from 
other nonleading Regge exchanges are much smaller, they are not completely neg- 
ligible and have been included in various analyses. 

The main problem, however, is posed by the sequence B w —> i — > PP in- 
volving inelastic rescatterring processes i — > PP. Arguments have been given 
that it is these inelastic processes that actually constitute the main source of soft 
FSI phases |j], § Q, f|. It has been also pointed out that nonzero inelastic- 
ity strongly affects the extraction of FSI phases in models based on quasi-elastic 
rescatterring. Thus, inelastic events affect model predictions even if rescatterring 
is of quasi-elastic type only. 

On the other hand, FSI phases are often attributed directly to short- distance 
(SD) quark-line diagrams in the hope that this will take into account all inelastic 
production phenomena. This belief persists despite justified scepticism about the 
dominance of short-distance QCD in FSI of B decays (see eg. |§). In fact, it 
is known that the resulting prescription strongly violates such tenets of strong 
interactions as isospin symmetry || (see also 0). The origin of the problem 
pointed out in ref. || is the lack of any correlation between the spectator quark 
and the products of b quark decay. By its very nature such correlation cannot be 
provided by SD dynamics. What must be involved here is a long-distance (LD) 
mechanism which ensures that quarks "know" about each other. Thus we are led 
to hadrons, hadron-level dynamics, and inelastic rescatterring effects. 

In this paper we perform a simplified analysis of corrections which should 
be introduced by inelastic rescatterring into the SD-based description of some 



nonleptonic decays of B mesons. Of course, any such analysis must be half- 
qualitative in nature, because - in the presence of many decay channels - it is 
well beyond our ability to take them accurately into account. Since it is believed 
that B decays into tctt and KK may provide some handle on the determination 
of angle a of the unitarity triangle, we shall concentrate on these decays: it 
is important to know the effects of FSI here. As shown in ref. ||, inclusion 
of coupled-channel quasi-elastic effects {nit — > KK) generates an effective long- 
distance penguin amplitude comparable in size to the short-distance one. One 
may expect that inelastic channels will also contribute to this effect. 

We start with an SD-based model of nonleptonic B decays. On the basis of 
standard tree-dominated mechanism for these decays, enriched with the related 
and well-established models of semileptonic B decays, we qualitatively estimate 
the types and the number of states produced in the first stage of the nonleptonic 
decay. As in other existing models, we take these states as composed of two (flying 
apart) resonances M\M 2 (Section 2). These resonances are assumed to rescatter 
into PP through Regge exchange. 

In order to provide the basis for an estimate of this rescatterring, we recall 
how in a Regge picture the unitarity relation involving M\M 2 — » PP and other 
% — > PP processes looks like. This enables us to make a rough estimate as to what 
part of all inelastic i —>■ PP processes is due to the M\M 2 — > PP transitions and, 
consequently, how much the situation deviates from the case of (quasi-)elastic 
rescatterring. Using a rough estimate for the contribution \{M\M 2 \S\PP)\ 2 from 
an average single inelastic intermediate channel M1M2, we estimate the number 
of inelastic channels involved (Section 3). 

In Section 4 we analyse the behaviour of the B — ► titt and B — ► KK amplitudes 
as a function of the number of intermediate states considered. We show that 
although effects due to each single intermediate state are small, the combined effect 
of all intermediate states is large. In particular, amplitudes of decays into KK 
become significantly larger than those estimated through short- distance penguin 
amplitudes, to the point of being comparable to the B — > nit amplitudes. 

Our conclusions are given in Section 5. 



2 B decays without FSI 



In short-distance approaches to nonleptonic weak B decays, the relevant ampli- 
tudes are usually expressed as sums of amplitudes corresponding to different types 
of quark diagrams (T - tree, C - colour-suppressed, E - ^-exchange, P - penguin, 
A - annihilation, PA - penguin annihilation). In this paper we concentrate on 
AS = AC = decays of B mesons which are initiated by a b —>■ uud transition. 
In these decays one expects that the dominant contribution comes from the tree 
diagram T, with the main corrections to it provided by the colour-suppressed 
and penguin diagrams C and P |J. Accordingly, neglecting contributions from 
other diagrams, the SD amplitudes ((P1P2) i\w\B°) for B° decays into a pair of 
octet pseudoscalar mesons P\P<i with total isospin / are expressed in the SU(3) 
symmetry case as 
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Inclusion of tti]s and r^s info our considerations is mandatory if we want to 
maintain SU(3) symmetry ||10|| . 



From the above formulas one may find SD amplitudes w^/ from state \B°) 
into states (R, I\ of given isospin / belonging to definite representations R of 
SU(3): 

wr,/ = O/W/ (2) 



with 



Wo 



Wi 



w 



T 



((tttt) 2 \w\B°) 

[((in^M^UAshM^ )] 

l((nn) \w\B°), ((KK) \w\B°), (( M8 ) |Hfi°}] 



(3) 



and the matrices O/ given by: 
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where rows correspond (from top to bottom) to 27 for O2; 8,27 for Oi; and 
1,8, 27 for O . 

One expects that SD amplitudes C and P constitute a 10-20% correction M. 
Indeed, for the factor r in the relation C = T/(3r), the short-distance QCD 
corrections give the value of r = {c\ + C2/3)/(3c2 +c{) ~ —3 (c\ ~ 1.1, c 2 ~ —0.25 
are Wilson coefficients), while estimates in ref. [11] yield \P/T\ in the range of 
0.04-0.20. 

Dominance of T amplitudes is expected to hold for other decays initiated by 
b — > uud as well. For larger invariant masses of du and ud systems these quark- 



level states should hadronize mainly as rm states [jT2| , |i3f . It is the rescatterring 
from these states to the final PP state that is of interest to us. Since we want 
to estimate rescatterring at hadron level, we need to know what hadronic states 
are produced in the first stage of the decay. There are two groups of hadrons 
produced: one comes from the decay W~ — > du, the other one originates from the 
recombination of the u quark with spectator d (Fig. 1). The values of invariant 
mass squares q 2 = s\ and s 2 (see Fig. 1 for definitions) are not large. The 
well-known probability distribution of 6-quark decay is 



v(q 2 ) = 2(1 -q 2 /m 2 b Y(l + 2q 2 /m 2 b ) 



(7) 



which falls for increasing q 2 , the average s\ being s7 m 7 GeV 2 . The values of s 2 are 
smaller. Estimates obtained in various models for semileptonic and nonleptonic 
decays yield I2 around 1.5 GeV 2 or so, with the distribution of ^Js2 extending 
to around 2.0 or 2.5 GeV |I2" , |I4"| , |T5| , [T6| , |TTt |j~8f . When such a low-mass quark- 
antiquark state hadronizes, a resonance is produced. In the ISGW2 model |16| for 



semileptonic decays, exclusive partial widths for the production of lowest-in-mass 
resonances have been predicted. The resonances considered were the ground-state 
mesons (pseudoscalars P: 1 1 Sq and vector mesons V: l 3 Si), the P-wave mesons 
(tensor mesons T: 1 3 P2, axial mesons A: 1 3 P\ and B: l 1 Pi, and scalar mesons 
S: 1 3 Po)j an d the 2 1 Sq and 2 3 S*i states. The highest (nonstrange) resonance mass 



explicitly considered in |Tj| is below 1.5 GeV. The total partial width for the 
production of all these resonances is 5-6 times larger than the partial width for 
the production of a pseudoscalar meson. Thus, the average partial width into a 
resonance is smaller by a factor of 0.7 or so than that for the production of a 
pseudoscalar meson. The resonances explicitly considered in [16[ do not saturate 
the inclusive decay rate which is still about two times larger |Tj|. Thus, several 
other resonances of masses below 2.0 or 2.5 GeV should be added to the list given 
above. Assuming that the average contribution from each one is similar to that 
just estimated, one expects that the number of types of "average" resonances 
produced should be of the order of 15. 

Similar or even larger number of resonance types is expected from the hadroniza- 
tion of the du created from the VT-boson. Thus, in nonleptonic B decays, apart 
from the PP state, many other resonance pairs must be produced: MiM 2 = VV, 
...,PA, PB, VA, VB, VS,...etc. Exact counting of the number of all these two- 
resonance states is not important for our purposes. However, it is fairly easy to 
give an estimate: limiting oneself to resonances of mass smaller than 2 GeV, this 
number will definitely be greater than 10, probably of the order of a few tens. 

Clearly, there is no hope that one can reliably calculate the contribution from 
rescatterring into PP from each of these intermediate states. However, one may 
try to estimate their overall contribution in an average way. For the first stage 
of the decay process, we will assume that the only important amplitude is the 
tree amplitude T and that this amplitude is approximately the same for all inter- 
mediate states considered (for low values of S\ that we shall be concerned with 
later, this is is indeed the case in Eq. ([/])). Although this may seem a very rough 
assumption for S2, we shall further see that our general results should be fairly 
independent of it as long as there is a rather large number of two-resonance states 
with production amplitudes scatterred around the average. The next question is 



how to describe transitions M\M 2 — > PP (or vice versa) in an average way. This 
is what we shall discuss in Section 3. 



3 Elastic scatterring and multiparticle produc- 
tion processes in PP interactions 

Elastic scatterring and multiparticle production processes are related to each other 
through unitarity of the S-matrix. Since B has spin J = 0, we shall work with 
Sj =0 sector only. In the following we shall suppress the subscript J and its value. 
In refs. || [H| it was shown that in the SU(3) symmetry case, with the effects 
of coupled channels included, one should work with states (PP)r, belonging to 
definite representations R = 1,8,27 of SU(3). The unitarity relation for PP 
scattering in the I = partial wave and in SU(3) representation R is: 

\((PP)n\S\(PP)n)\ 2 + T,\(( PP )n\S\k R )\ 2 = l (8) 

where, for given R, kn labels states different from (PP)n. Matrix elements oc- 
curing in Eq. fl8|) may be expressed in terms of Argand amplitudes a as follows 

((PP)n\S\(PP) R ) = l + 2ia((PP) R ) 

((PP)n\S\k R ) = 2ia(k R ) (9) 

Apart from the Pomeron, there are other Regge trajectories, whose exchange in 
the t-channel contributes to quasi-elastic scatterring (PP)-r —>■ (PP)r. When 
the leading non-Pomeron exchange-degenerate Regge trajectories p, / 2 , to, and a 2 
and their SU(3)-symmetric partners are taken into account, one obtains the / = 
partial wave amplitudes ||: 

a(( p P)27) = J-fa+^n 
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a{iPP)l) = iu + *tm±m?\ (10) 
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where for the leading non-Pomeron Regge trajectory we use 

a(t) = a(0) + o , t«0.5 + t, (12) 

ie. we have put a' = 1 GeV -2 and, consequently, in Eqs (|T0|) and further on both 
s and £ are in GeV 2 . 

The amplitudes a((PP)n) are independent of isospin / as shown in ||, where 
their sizes at s = m? B have also been estimated. The SU(3)-symmetric Regge 
residue R is fixed from experiment as || 

R/a' = -13.1mb GeV 2 = -33.6 (13) 

while for the Pomeron one has || 

P = 3.6 mb Gev 2 = 9.25. (14) 

Using the above values in Eqs ( |T0| ) one finds 

a((PP) 27 ) = -0.076 + 0.184i 
a((PP) 8 ) = +0.019 + 0.217z 
adPP)!) = -0.076 + 0.291i (15) 

while the leading non-Pomeron Regge contributions alone are: 

a((PP) 27 ,Reg) = -0.076 
a((PP) 8 ,Reg) = +0.019 + 0.033i 
a((PP) u Reg) = -0.076 + 0.107z (16) 

Contribution from elastic scatterring (the Pomeron in Regge language) is inde- 
pendent of R: 

adPP)^ Pom) = -!— P = 0.184 i (17) 

167T 



(cf. a = 0.17 % in ref. 0). Omitting the Reggeon-Pomeron interference term, the 
value of the contribution from the leading non-Pomeron Regge exchange to the 
unitarity relation of Eq. (|J), after averaging over representations R, is equal to 

\{PP\S Reg \PP}\ 2 = A\a(PP; Reg)\ 2 « 4|0.08| 2 = 0.025 (18) 

Neglecting Reggeon contribution to PP — > PP one obtains from Eq. (Q): 

(l-£) + y £\(PP\S\k)\ a = l (19) 

One may conjecture that contributions to the above sum from Reggeon-exchange- 
induced processes PP —>■ M\M 2 (where Mi denotes low-lying resonance) will be 
of size similar or smaller than \(PP\SR eg \PP)\ 2 . Thus, if all inelastic channels k 
were two-resonance states MiM 2 , one would obtain from Eq. (|19"D the number of 

ntot = £ \(PP\S\M 1 M 2 )\ 2 /0.025 = I 1 - f 1 - £-) J /0.025 « 25 (20) 

as the number of states contributing to the unitarity relation. This should be 
compared with the estimate of a few tens obtained in the previous section for the 
number of two-resonance states produced in weak decays of B meson. Of course, 
the estimate of Eq. (|20D is probably too low: contributions from transitions 
PP — > M\M 2 for heavier resonances Mi are likely to be smaller and the total 
number of states may be larger. For example, with average value of a going down 
by a factor of 0.6-0.7 from 0.08 to 0.05 (Section 2), the above estimate of the 
number of average quasi-two-body states increases from 25 to over 60. On the 
other hand, two-resonance states need not saturate Eq. flT9|) . Thus, there are two 
important questions which should be answered: 

1. how many of states k in the unitarity relation of Eqs fl3|,|l9"|) are indeed of 
the form MiM 2 , so that they can contribute to rescatterring in B decays, 
and 

2. can the effect of all of these states be described in some average way. 

At present, elastic and quasi-elastic contributions from long-distance FSI in 
B decays are usually evaluated using the old language of Regge theory. This 

9 



language was used in the past also for the description of resonance and multipar- 
ticle production processes that are of interest to us. In particular, the question 
of the build-up of elastic scatterring (Pomeron exchange) as a shadow of inelas- 
tic multiparticle production processes (and thus the very content of the unitarity 
relation of Eqs (|8"|JK)|) ) was discussed extensively. Therefore, we must recall the 
essential elements of an approach which dealt with that problem. The approach, 
predominantly occupied with the issue of the unitarization of dual models (for 
reviews see |19|]), was based on general properties of dual string models and on 
phenomenological analysis of resonance and multiparticle production data. 

In this approach, multiparticle production processes occurring in hadron-hadron 
collisions at high energies are pictured as proceeding through the production of res- 
onances or clusters in a multiperipheral model (Fig. 2) with leading non-Pomeron 
Reggeons exchanged in between the clusters. It is the sum over all types and 
numbers of resonances produced in this way that saturates the S-matrix unitarity 
relation Eq. (^). For s above the inelastic threshold but still small enough, one 
may limit oneself to the production of just one pair of resonances M1M2. With in- 
creasing energy s, the number of resonances produced in a single collision increases 
on average. Although cross-section for a production of a particular number of res- 
onances goes down at sufficiently high energy, the sum of cross sections over all 
possible numbers of resonances remains approximately constant. This constancy 
of the total inelastic cross section is ensured by a sufficiently fast increase in the 
number of all possible quark-line diagrams, ie. in the number of all possible states 
k (and ways in which they are produced) 



In our case, at s = m 2 B = 28 GeV 2 the model [21, 2l| predicts that states 
composed of just two resonances are produced in the fraction of /2M ~ 50% 
cases approximately. Further 35% comes from production of three resonances, 
etc. Although these numbers are obtained in |2TJ for resonances of any mass, 
the main contribution comes from production of resonances of invariant mass 
squared smaller than 6 GeV 2 . (Contribution from the production of objects of 
mass rriM is suppressed as (m 2 ^) 2 ^ ^ 1 ) for larger tum p3j ). The average mass 
of a resonance produced may be estimated in various ways to be around Wlm = 
1.5 (1.7 GeV in ref. [p3|), in good agreement with mass rnu ~ 1.3 or so |l 
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expected for average yfs^ in the SD-based models of weak decays. Contributions 
from rescatterring of states with larger values of s\ will be suppressed because 
production of such states in PP collisions is not likely. Thus, in a rescatterring 
process k — > PP one may expect that the dominant contribution will indeed 
come from the rescatterring of states composed of two low-mass resonances M1M2. 
Translating the above expectation of a 50% share of M\M 2 states in the unitarity 
relation of Eqs (|§,|T^) into a number n 2 M of contributing channels which may 
connect to the state originally produced by the SD dynamics (ie. n 2M = f^M^tot) > 
we conclude that this number should be around n 2 M = 50% x 25 ~ 12 for average 
\a(M 1 M 2 )\ « 0.08 or n 2M = 50% x 60 « 30 for average \a{M x M 2 )\ « 0.05. The 
latter estimate is probably more realistic since the average size of a contribution 
from a single M\M 2 channel should diminish with growing resonance masses. 



4 B decays with FSI 

If one accepts that final state interactions cannot modify the probability of the 
original SD weak decay, it follows that vector W representing the FSI-corrected 
amplitudes is related to vector w of the original SD amplitudes through 

W = S 1/2 w (21) 

Indeed, in the basis of S- matrix eigenstates |A) the above equation reduces to 
W\ = e lSx w\, ie. the condition of unchanged probability (\W\\ = \w\\) admits 
Watson phases only. 

The S-matrix may be written in terms of the matrix A of amplitudes a: 

S = 1 + 2iA (22) 

We assume that we may treat the FSI-induced corrections to the SD decay am- 
plitudes in a perturbative fashion. This is in agreement with the ideas of the 
dominance of SD dynamics. If this assumption is incorrect, obtaining even half- 
quantitative predictions will be almost impossible (cf. ref. [0]). Although this 
assumption may be questioned, it has an important advantage: one may study 
what happens when the number of contributing two-resonance intermediate states 
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is increased to its expected share (/ 2 m = 50%). In agreement with the assumption 
of a perturbative treatment of rescatterring (ie. small contribution from A) , we 
expand the square root in S 1//2 = (1 + 2iA) 1 ^ 2 and keep only the first term. This 

leads to 

1 + S 
W^(l + iA)w= w (23) 

This is in fact the K-matrix prescription for the estimate of rescatterring effects 
0. In this prescription, final state interactions modify probabilities of the SD 
amplitudes as can be seen explicitly from Eq. ( p3|) written in the basis of S- 
matrix eigenstates. Our calculations will be based on Eq. (p3|) . 



In previous sections we considered a simplified picture of contribution from 
resonance pairs in which all amplitudes a(MiM 2 ) were equal in absolute magni- 
tude while their phases were arbitrary. Indeed, when one row of the unitarity 
condition (ie. Eq. (||)) is discussed, no knowledge of phases is needed. For the 
purpose of studies of CP- violation the question of phases is important, however. 
Therefore, we have to make a very rough estimate of the FSI phases appearing in 
the strong rescatterring amplitudes ((PP) R \A\(M 1 M 2 ) R ) in 

((PP)n\W\B) = ((PP) R \w\B) + t J2 ((PP)r\A\(M 1 M 2 ) r )((M 1 M 2 ) u \w\B). 

MiM 2 

(24) 
Thanks to CP invariance of strong interactions, the ((PP)-r\A\(MiM 2 )-r) ampli- 
tudes are symmetric (as is the S- matrix). In order to estimate them we have to 
recall what are the predictions of dual string models for the production of two 
resonances in high energy PP collisions. In the Appendix of ref. |22] it is shown 



that the dual string model predicts that the amplitude for the PP —>■ M\M 2 pro- 
duction through an uncrossed diagram of Fig. 3a will pick up a rotating Regge 
phase resulting from the expression 

{-s/{s lS2 )r^. (25) 

Similarly, for the crossed diagrams of Fig. 3b one has to remove the " — " sign in 
the above expression, ie. the amplitude is real. Thus, the phase-generating factor 
differs from the familiar one in PP — > PP scatterring (ie. (— s) a ^) only by a 
different scaling factor (sis 2 ) in the denominator. Such a dependence of Regge 
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amplitudes on the masses of produced resonances has been confirmed in analyses 



of experimental data [2|], [2 



At this point we have to take into account the fact that Regge amplitudes 
describe scatterring of two colliding resonances Mi and M 2 in a state of definite 
momenta into a similar PP state. In particular, the produced PP state is a super- 
position of partial waves, while we are interested in the S j = q sector of the S'-matrix 
only. Restriction to the J = sector is achieved by integrating the rescatterring 
amplitudes a{M\M2) with P;=o(cos#) (with I being the angular momentum of the 
PP pair) over the allowed range of scatterring angle 9, or over the corresponding 
range of momentum transfer t G (t max ,t m i n ). Angular momentum conservation 
will then admit states M\M 2 with total angular momentum J = only. This 
will have no effect on the assumed form of the SD decay amplitude since in the 
SD mechanism the decaying b quark does not "know" about the spectator quark, 
and, consequently, about the value of J for the whole system. 

One calculates that for s\,S2 <C s the minimal value of t is t m i n ~ — sis-z/s, 
while for t max one may assume t max ~ —00. Projecting the M\M 2 — > PP Regge 
amplitude onto the J = sector, ie. integrating Regge expressions over t G 
(tmax, tmin), we obtain, up to a common overall normalization factor, the following 
phase factors: 

1. for the uncrossed diagram (Fig. 3a) 

z a ^exp(- (Inz-i7c)/z-i7ca(0)) , . 

9inei(z) = 7 . -^ (26) 

(In z — m) 

2. for the crossed diagram (Fig. 3b) 

*°(°)exp(-ln(*)/*) 

hnelW- ln( ^ [Zf) 



where z = sl(s\s<i). 



The overall normalization of contributions from amplitudes M\M 2 — > PP was 
fixed in Section 3. 

For the Regge description to be valid, the value of s/(sis 2 ) must be large. 
With s = m 2 B ~ 28 GeV 2 , the product S1S2 should not be greater than (s\S2)max ~ 
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6 GeV 4 , possibly 9 GeV 4 , corresponding to the minimum value of s/(sis 2 ) equal 
to 4 ±0.8. For s\ equal to s 2 , this corresponds to the maximum value of resonance 
mass being around 1.6-1.8 GeV. These are still reasonable numbers when com- 
pared with our previous estimates of Tom m 1.5. In Fig. 4 we show dependence 
of the phase of gi ne i{Tn%/(siS2)) on Sis 2 for 1 < S\S 2 < 9. If one approximates Si 
by their average values of around (1.5 GeV) 2 , one finds that the average value of 
the product Sis 2 is close to 5 GeV 4 ". From Fig. 4 we see that for Sis 2 ~ 4.7 GeV 4 
the phase of <ft ne / is zero. Although the phase of gi ne i at smaller and larger values 
of S1S2 deviates from zero quite significantly, these deviations are on average of 
the order of 20° and, consequently, it is still meaningful to talk about coherent 
superposition (ie. with approximately similar phases) of Regge contributions from 
various intermediate states. 

Because the quark-line structure of PP — > PP and M\M 2 — > PP amplitudes is 
the same, the eigenvalues with which the above phase factors enter into expressions 
for (MiM 2 )-r — ► (PP)n in definite SU(3) representations R are the same as those 
in amplitudes (PP)r —>■ (PP)-r. The latter were calculated in ref. ||. 

On the basis of ref. M we have therefore 

5(27, M 1 M 2 ) = a,rg{+2f ind {s/{ Sl s 2 ))) 
8(8,M 1 M 2 ) = 8£g(-~f inel (s/{s 1 s 2 )) + ~g ine i(s/(sis 2 ))) 

5(l,MiM 2 ) = Big(--f ine i(s/(s 1 s 2 )) + —g ine i(s/(s 1 s 2 ))j (28) 

The FSI-corrected amplitudes for P-meson decays into states in definite represen- 
tations R of SU(3) and with isospin I are 

W R)/ = (l+io((PP) R )w R>J +z J2 e i5(R ' MlM2) |a((MiM 2 ) R )|w R)7 (M 1 M 2 ) (29) 

M1M2 

where w R /(M 1 M 2 ) are SD weak decay amplitudes into M\M 2l all assumed ap- 
proximately equal to w R j (Eq. (g), Section 2). As in Section 3, we assume that 
all amplitudes |o((MiM 2 ) R )| are equal to some average value a. 

Although the phase of gi ne i changes over the range of corresponding Sis 2 , it 
is very instructive first to approximate it everywhere by a constant, namely its 
average evaluated at, say, 'sis 2 ~ (jTim) 2 ~ i.TGeV 4 , where gi ne i is real. The 
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approximate reality of average g ine i is a consequence of the particular value of s 
(being here equal to m 2 B ~ 28 GeV 2 ) and not an s-independent feature of the 
approach. At S1S2 = 4.7 GeV 2 we get 



5(1) « 0° 



5(8) « ±180° 



5(27) « 0° (30) 

The sum over all n%M two-resonance states will then yield the contribution from 
inelastic rescatterring: 



W K>I (inel) = i n 2M e i5 ^ a w Ri/ (31) 

Note that for fixed f 2 M the average amplitude 



1 " (l " fif ) hu 



\ 4n 2 



(32) 



is inversely proportional to ■Jni^M- Thus, contribution from inelastic events in 
Eq. (|31~D is proportional to yjniu'- the smaller the value of a, the larger the 
summed contribution from all two-resonance states, provided their contribution in 
the unitarity relation (Eq. (§)) is fixed by the same value of j 2 M- Obviously, this is 
a general feature of any perturbative treatment of rescatterring contribution from 
several channels, the reason being the linear nature of perturbatively treated FSI 
as compared to the quadratic nature of the unitarity relation. Thus, estimates of 
FSI effects given here are most likely estimates from below: although amplitudes 
a(MiM-2) corresponding to rescatterring from states composed of resonances of 
larger masses are expected to be smaller, their combined rescatterring effect should 
be relatively larger. 

The amplitudes for B decays into (tck)i, (KK)j are calculated from the inverse 
ofEq. (@): 

W/ = Of W Ri/ (33) 

In Table 1 we give the values of ((PP)i\W\B°) calculated from the above equation 
with the approximation of average phase (Eq. (|3l|) ) for PP = tttt, KK and 
n 2M = 12, 30 (a = 0.08,0.05 respectively). Predictions of the model without 
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inelastic rescatterring, ie. with quasi-elastic FSI (PP intermediate states) only, are 
also given for comparison. Amplitudes are given in units of input tree amplitude 
T. 

From Table 1 one can see that quasi-elastic FSI do not affect the decays 
B — ► (mr)i strongly. The amplitude moduli are somewhat smaller than those 
without the FSI. This is due mainly to the 1 + i [iP/(8ir)] factor originating 
from Pomeron exchange. Phase changes are small. For the B — > (KK) decays, 
quasi-elastic FSI affect the amplitude significantly: the amplitude driven in SD- 
dynamics by the penguin diagram (here vanishing) receives contribution from the 
coupled-channel-generated long-distance penguin (hereafter denoted Pld) @- Us- 
ing \((KK) \W\B°}/((7i7i) \W\B°}\ « y/3P LD /(2T) to estimate the effective LD 
penguin, the size of Pld is (as in ref. |§) of the order of 5% of the tree amplitude 
T, thus permitting significant interference effects with the short-distance penguin 
amplitude when the latter is taken from standard SD estimates. 

The main results of this work are given in the last two columns of Table 1. 

• One can see that the B — > (7T7r)j amplitudes, when compared with their 
estimates taking into account quasi-elastic FSI only, increase in absolute 
magnitude by a factor between 1 and 2. This is due to the additional 
contribution coming from rescatterring chain B — > M\M.2 — > im. 

• An important change can be seen in phase sizes: they are now one order 
of magnitude larger than in the quasi-elastic case. The origin of this effect 
is as follows. The W amplitude is composed of two parts: one (approx- 
imately real) is the SD amplitude w weakly suppressed by elastic (quasi- 
elastic) rescatterring, while the other contains the contribution from the 
inelastic MiM 2 — > 7r7r rescatterring. The latter part, being proportional to 
ia(MiM 2 )w, is mainly imaginary (for approximately real a). With a large 
contribution from inelastic rescatterring, the resulting phase must therefore 
be large. 

• Finally, we see that the B — > (KK) amplitude becomes much larger than 
in the case of quasi-elastic rescatterring only, and is dominantly imaginary. 
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Since the B — > (KK) amplitude is fed from no-hidden-strangeness states 
(MiM 2 )i=o, one might be tempted to compare this with ref. |]]. Namely, 
it was shown there for the case of a simple two-channel S-matrix that the 
phase of amplitude in channel 1 is large if the particle originally decays to 
channel 2. In our case, although a similar result holds, it is not general - 
it depends on the value of s at which amplitudes a(M 1 M 2 ) are evaluated. 
Note that inelastic rescatterring renders the B — > (KK)q amplitude larger 
than the SD penguin estimate: in fact, it is comparable to the B — > (7T7t)o 
amplitude. 

In our approach the full (FSI-corrected) amplitudes for B —* (KK) 1 decays 
vanish since: 

l)we have assumed that only the tree amplitude T is nonzero, and 
2) the rescatterring from the isospin \I = 1, I 3 = 0) state of (MiM 2 )i=i, which in 
principle might feed the final {KK)\ channel, is zero []. 

The dependence of the LD rescatterring-induced effective penguin on the isospin 
channel - with vanishing (large) effects in 1=1 (0) channel - should be compared 
with the SD mechanism which assigns the same size and phase to penguin ampli- 
tude P in both the (KK)o and [KK)\ decay channels. This is a general feature 
of long-distance dynamics: the size and phase of a quark-level diagram depend on 



what isospin (SU(3)) amplitude it contributes to [|T0 



Although the above expectation of FSI effect larger than that naively expected 
is fairly general, we have to analyse in some detail our assumption of replacing 
the Sis 2 -dependent phase of g ine i with an average (and vanishing) phase. Such 
an assumption would be well justified if a large fraction of M\M 2 states led to 
phases close to the average. Since for a given value of the product s\s 2 the 
phase is fixed, we need to know the density of two-resonance states as a function 



1 Vanishing of the rescatterring contribution can be seen as follows. The 1 = 1 state of M1M2 
is antisymmetric (ie. of the form (M + M~ — M~M + )/y/2), while the rescatterring contribution 
due to the uncrossed diagram of Fig. 3a is zero when evaluated in between the antisymmetric 
state I (Mi 1^2)1=1) and the symmetric state \(KK)i) (for the definitions of states, see Eq. (1) 
of ref. S). Rescatterring through diagram of Fig. 3b cannot change the type of quarks and, 
consequently, cannot induce a transition from the no- hidden-strangeness state \M1M2) into 
\KK). 
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of Sis 2 . Regge models |£| predict that the dependence of \a(MiM 2 )\ 2 on m 2 M . 
is proportional to l/m 2 M . for larger values of mjw r This fall of distribution for 
larger mass values should be clearly visible already at the beginning of the region 
where using Regge description becomes sensible, ie. at m 2 M . p» 4 GeV 2 or so. 
The distribution of s 2 in SD models vanishes even faster (it is fairly negligible 
above s 2 ~ (2.0 — 2.5 GeV) 2 ). Since we want to estimate corrections to SD 
models, direct use of mass distributions generated in SD-models might seem to 
be the simplest and most natural choice. Furthermore, as rescatterring of states 
with larger values of s\ is suppressed by the small size of the contribution from 
transitions M\M 2 — > PP, one might use the SD distribution of s 2 as that of Si as 
well. The problem is, however, that a significant part of the SD distribution of s 2 
(or si) corresponds to s 2 (or s\) < 1 GeV 2 . On the other hand, Regge amplitudes 
assume the form of (s/(sis 2 )) a( -^ only for Si,s 2 > 1 Gev 2 . For suy < 1 GeV 2 , 
one should replace si(2) with («') _1 = 1 GeV 2 . Thus, if Regge phases are to 
be reasonably evaluated we must use an appropriately modified distribution of 
s 2 . We model this situation in the simplest possible way: by assuming that the 
distribution of TriM k — \fsii vanishes below 1 GeV and above 2.25 GeV, while in 
between these values it is given by 

p{m M ) = 2.88 - 1.28m M (34) 

This yields the average value of "effective" vau equal to 1.42, which is also the 
value obtained from the ACCMM distribution of ^fs^ (see eg. ref. fll8l ) if con- 



tributions from rtiM < 1 GeV are replaced with contributions at tum — 1 GeV. 
Using p(rriM) as in Eq. (0), it is straightforward to evaluate the distribution p 2 
of Sis 2 . This distribution (Fig. 5) is peaked at Sis 2 ~ 2.2 GeV 4 , its median being 
at 3.6 GeV 4 , and the average sjsi at 4.4 GeV 4 . The tail above 9 GeV 4 contributes 
a few percent only. Thus, the assumption of sis^ ~ 4.7 GeV 4 used in our previous 
discussion appears quite reasonable. 

With two-resonance states k spread in s\s 2 according to distribution p 2 (larger 
values of k — l,..,n 2 M correspond to states with appropriately larger values of 
Sis 2 ), it is simple to analyse the predictions of the model numerically. We are 
interested in the question how the FSI effects change when heavier and heavier 
intermediate states k are included. For decays into states in definite representa- 
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tions of SU(3), the amplitudes of interest to us are therefore: 

W R)/ (n) = (1 + ia((PP) K )w n ,i + i £ e l5 ^ s ^ a w Rj/ (35) 

fc=i 

where n is the number of intermediate inelastic states considered (0 < n < ti^m) 
and S(H, S1S2) are given in Eq. (|28|) . 



In Fig. 6 we present predictions of the n-dependent version of Eq. ( p3|) for the 
absolute values and phases of the amplitudes of B decays into (tttt)o and (KK) 
states. These predictions are given as a function of the number n of intermediate 
states considered. We show the case with u^m — 12 (for n 2 M — 30 one obtains 
very similar plots with features discussed below being even more pronounced). 
The point most to the right in each plot (ie. at n — U2m) corresponds to a 
large value of S1S2, where Regge approximation breaks down. Consequently, this 
point should be discarded. One can see that for B — > (7T7r) the size of the 
amplitude does not depend very strongly on n and is close to the value of the 
input FSI-free amplitude |((7T7r)o|u>|-B°)|, which is 0.58. Furthermore, the FSI- 
induced phase is still relatively small (of the order of —5°). This cannot be said 
of the B — » (KK) process. Here, the absolute value of the amplitude grows 
fast with the increasing number of intermediate states. In addition, already at 
n = 4 or so, the phase becomes close to 100° , not far from the previous result 
of 93° (Table 1) obtained for constant phases. One can also see that when the 
number of intermediate states approaches its maximum allowed value, the size 
of the B — > (KK)o amplitude becomes significant when compared to the tree 
B — ► (7T7r)o amplitude. As the number of intermediate states increases, the long- 
distance-induced penguin amplitude starts to dominate over the SD one (which 
was estimated at 0.04 to 0.20 of the tree amplitude, |]lljj). Furthermore, assuming 
applicability of Eq. ([I]), one estimates from Fig. 6 that at large n the effective long- 
distance penguin amplitude Pld is around 0.6 to 0.8. This should be compared 
with the input or effective T amplitudes which are around 1.0. Thus, the long- 
distance effective penguin amplitudes become really large. 

Large rescatterring effects obtained above stem from an approximately co- 
herent superposition of contributions from several intermediate channels. One 
might argue that in the calculation of this paper sizes of individual rescatter- 
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ring amplitudes a (ie. contributions from individual intermediate channels) are 
overestimated. However, the trend of the results expected for the case of smaller 
amplitudes a may be seen from Table 1 and Eqs. ([H],[!2|) . The column of n 2 M = 30 
corresponds to a smaller value of average amplitude a = 0.05. The connection 
existing between the values of n^M to a stems from the assumption that at energy 
s = m 2 B a two-resonance state is produced in approximately half of all inelastic 
PP collisions. As long as this fraction (./W) is kept constant and the perturbative 
treatment of FSI is valid, smaller values of amplitudes a - after summing over all 
intermediate channels - result in rescatterring corrections to weak decays larger 
than naively expected. One may also try to estimate roughly the rescatterring 
contribution in a direct way (ie. without using the value of /W) by just adding 
rescatterring contributions with amplitudes a{M\Mz) assumed to be of the order 
of a(PP) or so (in agreement with experiment). Then, for the number of aver- 
age intermediate states taken as equal to the number of final states in standard 
models of SD decays, ie. of the order of 10 or more, one is bound to obtain a 
large rescatterring effect. Please note that our whole approach starts with the 
generally accepted features of the SD decay (and resonance production) mecha- 
nism and assumes that subsequent FSI may be treated perturbatively. We have 
shown that even in this case the corrections tend to become large. Of course, if 
they are too large, the whole perturbative scheme of their estimation (as well as 
the SD mechanism for the description of B decays) ceases to be viable. 



5 Conclusions 

We have discussed multichannel inelastic rescatterring effects in B decays into tctt 
and KK. Generally accepted features of short-distance decay mechanism have 
been assumed as part of our input. These assumptions included estimates of the 
types and number of resonances produced in two-resonance states initiated by the 
b — » uud transition. Rescatterring of these two-resonance states into the final 
state consisting of one pair PP of pseudoscalar mesons was evaluated under the 
assumption that such FSI may be treated perturbatively. The basis for this eval- 
uation was provided by existing knowledge about how the inelastic multiparticle 
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(resonance) production in PP collisions is correlated with elastic PP (Pomeron 
exchange) scatterring. This knowledge permitted us to estimate that at s = m 2 B 
PP scatter into a two-resonance M\M.2 state in approximately 50% of cases. Thus, 
the total size of rescatterring from M\M.2 to PP was fixed. 

Using the Regge model for the description of the M1M2 — > PP processes, we 
have shown that the rescatterring contributions from the individual intermediate 
channels add approximately coherently. As a result, the combined effect of rescat- 
terring through many two-resonance intermediate states was shown to be quite 
large. This was demonstrated under the assumption that FSI may be treated per- 
turbatively. If that assumption is overoptimistic, reliable estimate of (presumably 
even larger) FSI effects will almost certainly be much more difficult. 

In our calculations the amplitude for the decays B — > (KK) was induced 
by LD rescatterring from no-hidden-strangeness isospin states produced via 
short-distance tree diagram. This FSI-induced amplitude was shown to be larger 
than its short- distance penguin counterpart. The phase of the LD amplitude was 
estimated to be around 100°. 
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FIGURE AND TABLE CAPTIONS 

Fig.l 

Tree diagram (T) for B decay 

Fig.2 

Multiperipheral production of multi-resonance state \M 1 M 2 ...M V ) through Regge 

exchanges Rk 

Fig.3 

Quark-line diagrams for production of two- resonance state \MiM 2 ): (a) uncrossed 

Reggeon exchange, (b) crossed Reggeon exchange 

Fig.4 

Dependence of 5wK m s/( s i s 2)) phase on sis 2 

Fig.5 

Distribution of S1S2 in model defined in text 

Fig. 6 

Dependence of FSI-induced effects on the number of intermediate channels in- 
cluded. 



Table 1. 

Effects of inelastic rescatterring on ((PP)i\W\B°) amplitudes in average phase 

approximation. Amplitudes are in units of input tree amplitude T. 
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Table 1 



(PP)i amplitude 

modulus/phase 


noFSI 


quasi-elastic FSI 


inelast 
n 2M = 12 


icFSI 
n 2M = 30 


(7T7T) 2 \W\ 

&rg(W/w) 


0.41 
0° 


0.33 

-5.3° 


0.49 

47° 


0.67 
60° 


(vr7r)o \W\ 

&rg(W/w) 


0.58 
0° 


0.44 
-1.4° 


0.44 
-16° 


0.61 

-17° 


(KK) \W\ 
arg(W) 





0.028 
223° 


0.36 
93° 


0.57 
93° 
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